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Inspired by recent measurements with the CLAS detector at Jefferson Lab, we perform a self- 
consistent analysis of world data on the proton structure function g\ in the range 0.17 < Q 2 < 
30 (GeV/c) 2 . We compute for the first time low-order moments of gi and study their evolution 
from small to large values of Q 2 . The analysis includes the latest data on both the unpolarized 
inclusive cross sections and the ratio R = ol/pt from Jefferson Lab, as well as a new model for 
the transverse asymmetry A 2 in the resonance region. The contributions of both leading and higher 
twists are extracted, taking into account effects from radiative corrections beyond the next-to- 
leading order by means of soft-gluon resummation techniques. The leading twist is determined with 
remarkably good accuracy and is compared with the predictions obtained using various polarized 
parton distribution sets available in the literature. The contribution of higher twists to the g\ 
moments is found to be significantly larger than in the case of the unpolarized structure function 
F 2 . 
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I. INTRODUCTION 

One of the fundamental characterizations of nucleon 
structure is the distribution of the nucleon spin among its 
quark and gluon constituents. The classic tool for study¬ 
ing the quark spin distributions experimentally has been 
inclusive lepton scattering off polarized protons and neu¬ 
trons. These experiments have determined the (j\ struc¬ 
ture function of the nucleon, which, in the framework of 
the naive Quark-Parton Model (QPM), is proportional 
to the difference between the distributions of quarks with 
spins aligned and anti-aligned to the nucleon spin. Sur¬ 
prisingly, one finds that only 20-30% of the proton spin 
is carried by quarks - an observation which came to be 
known as the “proton spin crisis”. Considerable effort, 
both experimentally and theoretically, has subsequently 
gone into understanding where the remaining fraction of 
the proton spin resides - see Ref. i'l'j for recent reviews. 

In terms of kinematics, most of the experimental study 
has been focused on the high-Q 2 region, where the QPM 
description is most applicable, and in the region of in¬ 
termediate and small Bjorken-x, which is important for 
evaluating parton model sum rules such as the Bjorken 
sum rule. Qualitatively new information on the proton 
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spin structure can be obtained by studying the gi struc¬ 
ture function in the region of large Bjorken-x, at moder¬ 
ate values of the squared four-momentum transfer Q 2 , in 
the range from 1 to 5 (GeV/c) 2 . Such a kinematic region 
is characterized by the presence of nucleon resonances 
which contribute to higher-twist effects in the structure 
functions. 

According to the operator product expansion (OPE) 
in QCD, the Q-evolution of structure function moments 
can be described in terms of a 1 /Q 2 , or twist, expan¬ 
sion, where the leading twist [0(1) in 1 /Q 2 ] represents 
scattering from individual partons, while higher twists 
[0(1/Q 2 ) and higher] appear due to correlations among 
partons. The inclusion of the contribution from the nu¬ 
cleon resonance production regions is a relevant point of 
our study, because resonances and Deep Inelastic Scat¬ 
tering (DIS) are closely related by the phenomenon of 
local quark-hadron duality |2j, -3, 4J. The latter has been 
extensively investigated at Jefferson Lab (JLab) for the 
case of the unpolarized structure function F 2 of the pro¬ 
ton Hi In the polarized case, the contribution of the 
A(1232) resonance makes the analysis rather more inter¬ 
esting: since this resonance gives rise to a negative con¬ 
tribution to the g 1 structure function, while g\ at high 
Q 2 is positive, one expects a breaking of local duality to 
occur in the A region at least up to several (GeV/c) 2 0. 

In this paper we report the results of a self-consistent 
extraction of the proton structure function g\ (x, Q 2 ) and 
its moments from the world data on the longitudinal po- 
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larization asymmetry A \\. The extraction is based on a 
unique set of inputs for the structure function F 2 , the 
ratio R = and the transverse asymmetry A 2 . The 

complete data set measured at Jefferson Lab J| j§, fiol . 
which covers the entire resonance region with high pre¬ 
cision, allows for the first time the (Revolution of the 
g 1 moments to be accurately evaluated up to n = 7. 
The results for the first moment have been presented in 
Ref. El , where the twist-four matrix element was ex¬ 
tracted, and the proton’s color electric and magnetic po¬ 
larizabilities determined. Here we give the details of our 
analysis for all the moments up to n = 7. 

In Section HQ we describe the OPE framework of the 
moments analysis for the polarized structure function g\. 
In Section UTTI we discuss the extraction of g± from the lon¬ 
gitudinal asymmetry A\\. The evaluation of the moments 
of g 1 and their uncertainties is presented in Section cm 
and the extraction of both leading and higher twists is 
described in Section m Finally, conclusions from this 
study are summarized in Section fVl 


II. MOMENTS OF THE STRUCTURE 
FUNCTION gi 

The complete (Revolution of the structure functions 
can be obtained using the OPE G2 of the time-ordered 
product of the two currents which enter into the virtual 
photon-nucleon forward Compton scattering amplitude, 

T \J{Z) J(0)] = /nR 2 ) Z^Z^...Z^ 0“ 1M2 ... Mb (1) 

n,a 

where CU lAt2 ... Mn are symmetric traceless operators of di¬ 
mension and twist — n, with a labeling dif¬ 

ferent operators of spin n. In Eq. 0, /“(— z 2 ) are coeffi¬ 
cient functions, which are calculable in perturbative QCD 
(pQCD) at short light-cone distances z 2 = (ct) 2 — z 2 ~ 0. 
Since the imaginary part of the forward Compton scat¬ 
tering amplitude is simply the hadronic tensor contain¬ 
ing the structure functions measured in DIS experiments, 
Eq. 0 leads to the well-known twist expansion for the 
Cornwall-Norton (CN) moments of gi(x, Q 2 ) pin, 

M n CN (Q 2 )= [ dx x 71 - 1 g?(x,Q 2 ) 

Jo 

00 / 2 \ («- 2)/2 

= E E n 4^a s (Q 2 )] O nK (fi) (gj) (2) 

k=2,4 ••• ' 

for n = 1,3, 5,.... Here g is the renormalization scale, 
O nK (/r) are the (reduced) matrix elements of operators 
with definite spin n and twist k. containing information 
about the nonperturbative structure of the target, and 
E nK (gL : Q 2 ) are dimensionless coefficient functions, which 
can be expressed perturbatively as a power series of the 
running coupling constant a s (Q 2 ). 

In the Bjorken limit ((R v —> 00 , with x = Q 2 /2Mv 
fixed, where v is the energy transfer and M the nucleon 


mass), only operators with spin n contribute to the n- 
th CN moment ©. At finite Q 2 , however, operators 
with different spins can contribute. Consequently the 
1/Q 2 expansion of the CN moment _AR N ((R contain 
in addition target-mass terms, proportional to powers of 
M 2 /(R which are formally leading twist and of pure 
kinematical origin. It was shown by Nachtmann 1 1 fill in 
the unpolarized case, and subsequent ly g eneralized to the 
polarized structure functions in Ref. 1141 , that even when 
M 2 /Q 2 is nonzero, the moments can be redefined in such 
a way that only spin-n operators contribute to the n-th 
moment. This is achieved by defining the “Nachtmann 
moments” of g\ as 



x n 2 M 2 x 2 £ 
£ (n + 2) 2 Q 2 x 



- 92 (x,Q 2 ) 


M 2 x 2 4 n i 
Q 2 n + 2 J 


( 3 ) 


where £ = 2x/ ^1 + i/l + 4 M 2 x 2 /Q 2 ^ is the Nachtmann 
scaling variable. Note that the evaluation of the polar¬ 
ized moments M n (Q 2 ) requires the knowledge of both 
structure functions g\ and g 2 . In the DIS regime the 
contribution of go to Eq. © turns out to be typically 
small (see Ref. |7j). On the other hand, in the nucleon 
resonance production region the impact of g 2 is expected 
to be more significant, and here the lack of experimen¬ 
tal information on the structure function g 2 can lead to 
systematic uncertainties. 

Since the moments in Eq. are totally inclusive, the 
integral in the right hand side of Eq. 0 contains also 
the contribution from the elastic peak located at x = 1, 


,2, xf ^ ,„2 ,Ge{Q 2 ) + tGm{Q 2 ) 


gf(x, Q 2 ) = S(x — i)G m (Q 2 )- 


2(1 +t) 


^ » u, ^ 2 R b (Q 2 )-Gm(Q 2 ) 


g 2 {x,Q ) = 5{x - 1 )tG m (Q )- 


2(1 + r) 


( 4 ) 

( 5 ) 


with Ge ( Gm) the proton electric (magnetic) elastic form 
factor and r = Q 2 /AM 2 . 

Note that the structure function moments include the 
resonance production region at low Q 2 and high x, which 
would be otherwise problematic to include in a twist anal¬ 
ysis performed directly in x-space. In addition, since 
target-mass corrections are by definition subtracted from 
the moments 0 , the twist expansion of the Nachtmann 
moments M n (Q 2 ) directly reveals information on the 
nonperturbative correlations between partons, without 
relying on specific assumptions about the x-shape of the 
leading twist. 

For the leading twist contribution [k = 2 in Eq. 0 ], 
one finds the well-known logarithmic Q 2 evolution of 
both singlet and non-singlet moments. However, if one 
wants to extend the analysis to small Q 2 and large x, 
where the rest of the perturbative series becomes sig¬ 
nificant, some procedure for the summation of higher 
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orders of the pQCD expansion, such as infrared renor- 
malon models Jjdl E3 or soft-gluon resummation tech¬ 
niques El 19, l20l| . has to be applied. For higher twists, 
k > 2, the power-suppressed terms are related to quark- 
quark and quark-gluon correlations, as schematically il¬ 
lustrated in Fig.^ and should become important at small 

Q 2 - 


FIG. 1: Example of a twist-4 diagram appearing in the OPE 
of polarized structure function moments. Arrows indicate the 
spin projections of the particles. 

The evaluation of the Nachtmann moments 0 from 
available data in the range 0.17 < Q 2 < 30 (GeV/c) 2 
will be described in the next section. The OPE analy¬ 
sis of such experimental moments will allow us to extract 
simultaneously both the leading and the higher twist con¬ 
tributions. A precise evaluation would permit a compar¬ 
ison of the leading twist with the QCD predictions ob¬ 
tained from lattice simulations or with nonperturbative 
models of the nucleon. 


III. DATA ANALYSIS 

The data analysis was performed starting from mea¬ 
sured longitudinal proton asymmetries Ay, which were 
converted into the structure function g\ using consistent 
values of the ratio R = olI^t and the structure func¬ 
tion Fj, as well as of the transverse proton asymmetry 
Ai. Our procedure is described in detail in the following. 

A. Asymmetry Database 

All available world data on the longitudinal and 
transverse asymmetries, Ay and A i , w ere collected 
from Refs [I EH El 0 0 0 IS 0 0 and 
Refs. |2l| | .[22| p [23l | b.[2fi | b. respectively. The full data set 
of A || consists of two subsets corresponding to the reso¬ 
nance 0 and DIS regions 00000000. 
The kinematic coverage of the experimental data is 
shown in Figs. |3 and 0 for Ay and Aj_, respectively. It 
can be seen that the resonance region is completely cov¬ 
ered by the Ay data up to Q 2 = 2.5 (GeV/c) 2 with the 
inclusion of recent high quality data from CLAS [Ij.In 
contrast, the A± asymmetry is poorly determined in the 
resonance region. The lack of data on A± here becomes 
problematic because of the prominent role of the higher 
twist contributions at large values of x. 



x 


FIG. 2: Kinematics of Ay world data from Refs. 100 
000000 (different symbols indicate different 
experiments). 



X 


FIG. 3: Kinematics of A± world data from 

Refs. ’211,^21. l2iil b.l2fib (different symbols indicate dif¬ 
ferent experiments). 


B. Extraction of the Structure Function g\ 

In order to extract the structure function g\ from the 
data collected in our data base, one needs additional ex¬ 
perimental inputs for the structure function F \, the ratio 
R , and the transverse asymmetry A 2 . Indeed, the struc- 
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ture function gi is given by 
gi{x,Q 2 ) = 


2 . _ F 1 (x,Q 2 ) j A\\{x,Q 2 ) 


1 + 7 2 1 D 
+ (7 - v)A 2 (x,Q 2 ) 


( 6 ) 


with 

2Mx jToy r _ (7I 

v^’ E — eE’’ 1 + eR(x,Q 2 )' U> 

where E and E’ are the incident and scattered electron 
energies and e is the virtual photon polarization. The 
ratio R entering above was taken from the parameter¬ 
ization given in Ref. Ilfll for the resonance production 
region, while in the DIS domain the fit R1998 (2l| was 
used. 

Since the main goal of our analysis is a model inde¬ 
pendent extraction of the moments of gi, the structure 
function Fi(x,Q 2 ) has been obtained directly from ex¬ 
perimental data. This has been possible because of the 
large amount of high quality data on the inclusive elec¬ 
tron scattering cross section da/dildE 1 and on the struc¬ 
ture function F 2 , covering both the resonance and DIS 
regions (for the list of data used see Ref. H). Therefore, 
for each point of the measured longitudinal asymmetry 
A\\ we can find several nearby points with either F 2 or 
the inclusive cross section known from experiments. For 
the interpolation of iq (x. Q 2 ) points, a simple procedure 
has been used, which is described below. 

Having a data point with the measured A\\ at some 
fixed Xq and Q 2 , we search in the combined database on 
the inclusive cross section dcr/dQdE ' and the structure 
function F 2 for several nearby experimental points. The 
search procedure chooses a rectangular bin around the 
point with coordinates (xq, Q q ) of such a size that the 
selected area contains a number N of experimental points 
either from da/dfldE' or from F 2 . The procedure then 
selects only those configurations whose number of points 

A rn'i/n N <: ' Nrnaxi where Nmin — 2 and Nrnax — 6 in 

the resonance region and N m i n = 1 and N max = 4 in the 
DIS case. Once a number of configurations have been 
collected (no more than 20 sets), the procedure looks for 
a minimum in the sum of the path integrals from each 
point ( Xi,Q 2 ) of measured da/dfldE' or F 2 to the bin 
center ( x 0 , Q q ), 


S(xq, Q 0 ) 


1 

NF^Q 2 ) 



dl |.Fi(ie, Q 2 )\ (8) 


where the integral over dl is taken along a straight line 
connecting the point (xt,Q 2 ) to the bin center (xq,Qq). 
The structure function Fi(x,Q 2 ) in this integral is con¬ 
structed using the fits of F 2 from Ref. [H] and of R from 
Ref. Ell in DIS, while in the resonance production region 


F\ is taken directly from Ref. 0- The configuration se¬ 
lected is that which minimizes the function S(x,Q 2 ) in 
Eq. ©. 

From Fig. |3 and also from Fig. 1 of Ref. Q, one can see 
that in the resonance region, which is covered by the data 
from Ref. Q, the interpolation distances are very small, 
thanks to the measurements of inclusive cross section in 
the same kinematic range 11 A set of experimental 
points of dcr/dQdE' or F 2 identified above is converted 
to the structure function F\ according to 


Fi(x,Q 2 ) 


and 


Fi(x,Q 2 ) 


MQ 2 E 1 — e da 
2a 2 E' 1 + tR{x,Q 2 ) dVLdE' 


1 + 4AI 2 x 2 /Q 2 
2x(l + R(x,Q 2 )) 


F 2 (x, Q 2 ) . 


( 9 ) 


( 10 ) 


All the F\ points obtained within the given bin are aver¬ 
aged together with their aq and Q 2 coordinates, 


i 


Fi{xj,Q 2 ) 
S 2 Fl (xi,Q 2 ) ’ 


1 ^ x i 

52 i 



(11) 

( 12 ) 
(13) 


where 



(14) 


and S Fl is the statistical error of F\. The mean value of 
F\(x 1 Q 2 ) is then corrected by the bin centering correc¬ 
tion using the models of Refs. 0 HE!- The value of 
the correction turns out to be very small with respect to 
statistical and systematic errors of the A\\ data. Never¬ 
theless, the correction value has been propagated in the 
total systematic error obtained for F\. 

Once the transverse asymmetry A± is known, A 2 can 
be determined according to 


where 


A 2 


1 

(1 + ^/0 


C£[ 

D 


A± 

d 


(15) 


d = 



C = v 


l + e 
2e 


(16) 


Since there are no experimental data on A± in the reso¬ 
nance region (see Fig.|3), we consider several models: 


• The model-independent constraint provided by the 
Soffer limit l3lj : 


\A 2 \< ] J^±1r. (17) 

This inequality is exact and, provided A\ and 
R(x,Q 2 ) are measured, gives unambiguous limits. 
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• Since it was shown in previous experiments that Ai 
is in fact much smaller than the Softer limit (2^ , one 
can simply assume A -2 = 0, with possible deviations 
from zero included in the systematic error. 

• In the present analysis we use a somewhat more 
sophisticated model for A 2 which is described in 
detail in Appendix 1X1 

The Q 2 dependence of gi(x, Q 2 ) at x = 0.38 — 0.42 is 
shown in Fig. fusing different assumptions about A 2 and 
F\ , which provides an estimate of the systematic errors. 
The ranges and the averages for the various sources of 
systematic errors on g\ are collected in Table Q] 



FIG. 4: Q 2 dependence of the structure function t/i at 

x = 0.38 — 0.42 obtained from the data in Refs. [3 I 2 H 
M M El Ei El ® using the procedure described in the 
text. Open squares represent central values obtained with 
the A 2 model described in Appendix m while the filled tri¬ 
angles indicate upper and lower Softer limits. The upper 
hatched area represents the difference between g 1 data points 
extracted with two different parameterizations of R Mfl 1171 ; 
middle hatched area F^ IM shows the difference between g\ 
data points extracted using two different parameterizations 
of F\ If}, ^3; lower hatched area F\ >M shows the difference 
between g\ extracted from the F\ parameterization and data 
interpolation as described in the text. 


C. Moments of the Structure Function gi 

As discussed in the introduction, the final goal of our 
data analysis is the evaluation of the Nachtmann mo¬ 
ments of the structure function g\. The total Nachtmann 


TABLE I: Range and average of systematic errors on g\ (ab¬ 
solute value). 


Source of uncertainties 

Variation range 

Average 


10' 4 - 0.14 

0.015 

Fi 

10~ Y - 1.7 

0.014 

ol/ot 

10~ 4 - 0.015 

0.002 

A.2 

1(T Y - 0.015 

0.004 

Total 

10” 4 - 1.7 

0.025 


moments were computed as the sum of the elastic (M® 1 ) 
and inelastic ( M™) moments, 

M n {Q 2 ) = M®HQ 2 ) + AC(Q 2 ) . (18) 


The contribution from the elastic peak can be calculated 
by inserting Eqs. into Eq. 0 , 


m:\q 2 ) 


ir\1\ f Ge{Q 2 ) + tG M (Q 2 ) 
y G m (Q )j- — - 


n 2 M 2 , 2 
X (n + 2) 2 Q 2 U 

Gm(Q 2 )~G e (Q 2 ) n \ 

-m-T-TA8e 1 > 


(19) 


where £ el = 2/(1 + yj\ + 1/r). 

The evaluation of the inelastic moments M™ involves 
the computation at hxed Q 2 of an integral over x. In 
practice the integral over x was performed numerically 
using the standard trapezoidal method in the program 
TRAPER H3. 

The Q 2 -range from 0.17 to 30 (GeV/c) 2 was divided 
into 24 bins increasing logarithmically with Q 2 . Within 
each bin the world data were shifted to the central bin 
value Qq using the fit of gf(x, Q 2 ) from Ref. jJ, which 
covers both the resonance and DIS regions, 

9i{x,Qo) = gi(x, Q 2 ) + (gf(x,Q 2 ) -gf(x,Ql)) . (20) 

The difference between the actual and bin-centered data, 

Sgt n \x, Q 2 ) = \gf{x, Ql) - gf{x, Q 2 )\ , (21) 

is added to the systematic error of g 1 in the Nachtmann 
moments extraction procedure. As an example, Fig. 0 
shows the integrands I n (x,Q 2 ) of two of the low-order 
moments as a function of x at fixed Q 2 . The significance 
of the large-a; region for higher moments can be clearly 
seen. 

To obtain a data set dense in x, which reduces the error 
in the numerical integration, we performed an interpola¬ 
tion at each fixed Qq when two contiguous experimental 
data points differed by more than V. The value of V 
depends on kinematics: in the resonance regions, where 
the structure function exhibits strong variations, V has 
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FIG. 5: Integrands of the Nachtmann moments at Q 2 = 
1 GeV 2 for the n = 1 (upper) and the n = 3 (lower) moments. 


to be smaller than half of the resonance widths, and is 
parameterized as V = 0.03 M 2 IQ 2 . Above the reso¬ 
nances, where g\ is smooth, to account for the fact that 
the available x region decreases with decreasing Q 2 , we 
set V = 0.1. Finally, in the low x region (x < 0.03) where 
the g\ shape depends weakly on Q 2 , but strongly on x, 
we set V = 0.005. 


where Sf^^Xj, Qq) is the 9i statistical error and 


~\xi—x a \<A 

5 n (qd = y - ¥ 

L i (^te.QS)) 

\xj — Xb|<A “I —1/2 

+ E 1 , 

i (sy( Xj ,Q 2 )y 


(24) 


is the statistical uncertainty of the normalization. There¬ 
fore, the statistical error of the moments calculated ac¬ 
cording to the trapezoidal rule (3(1 was increased by 
adding the linearly correlated contribution from each in¬ 
terpolation interval as 


£n+l 

dx ~^~9f Qo) 

X x n2 M2x2 . (25) 

£ {n + 2) 2 Qq x 

Since we average the difference gi(xi, Qq) — gf(xi, Qq), 
A is not affected by the resonance structures, and its 
value is fixed to have more than two experimental points 
in most cases. Therefore, A is chosen to be equal to 0.15. 

To fill the gap between the last experimental point and 
one of the integration limits ( x a = 0 or Xb = 1) we per¬ 
formed an extrapolation at each fixed Q 2 using gf(x, Qq) 
including its uncertainty given in Ref. 7]. The results, 
together with their statistical and systematic errors, are 
presented in Table un 


C rm (Qo) = MQo) 


D. Systematic Errors of the Moments 


To fill the gap between two adjacent points x a and Xb, 
we used the interpolation function g[ nt (x, Qq), defined as 
the parameterization from Ref. @ offset to match the ex¬ 
perimental data on both edges of the interpolating range. 
Assuming that the shape of the fit is correct, one has 


g\ nt (x, Ql) = p(Ql) + gf (x, Ql) , (22) 


where the offset p(Qq) is defined as the weighted average, 
evaluated using all experimental points located within an 
interval A around x a or Xb'- 


p{Ql 


= Fn(Qo 

\xj—x b \<A 


'\xi—x a \<A 


51(24, Qo) - 9i (24, Qo) 


E 


y 

i {siT(xi,Q 2 0 )Y 

gi (xj , Qq) — gf (xj , Qq) 

(sy( Xj ,Q 2 )) 2 . 


(23) 


The systematic error consists of experimental uncer¬ 
tainties in the data given in Refs. H, |(01152, H2 |3> EE 
[HElEi and uncertainties in the evaluation procedure. 
To estimate the first type of error we have to account for 
using many data sets measured at different laboratories 
and with different detectors. In the present analysis we 
assume that different experiments are independent and 
therefore only systematic errors within a particulardata 
set are correlated. 

An upper limit for the contribution of the systematic 
error from each data set was thus evaluated as follows: 

• we first applied a simultaneous shift to all experi¬ 
mental points in the data set by an amount equal 
to their systematic error; 

• the inelastic n-th moment obtained using these dis¬ 
torted data MnY (Q 2 ) is then compared to the orig¬ 
inal moments M™(Q 2 ) evaluated with no system¬ 
atic shifts; 
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TABLE II: The inelastic Nachtmann moments for n = 1,3,5 and 7 evaluated in the interval 0.17 < Q 2 < 30 (GeV/c) 2 . The 
moments were evaluated for Q 2 bins with more than 50% data coverage. The data are reported together with the statistical 
and systematic errors; the low-a: extrapolation error is given for the first moment only (last number in the second column). 


Q* [(GeV/crj 

Mi(Q") x KT J 

M 3 (Q") x 10" 4 

M 5 (Q z ) x 10 -b 

M 7 (Q z ) x 10“° 

0.17 

-27.1 ± 7 ± 12 ± 6 

16.8 ± 2.5 ± 5 

-8.5 ± 1 ± 2.5 

-4.8 ± 0.6 ± 1.3 

0.20 

-23.0 ± 5 ± 9 ± 6 

-17.0 ± 2 ± 4 

-8.4 ± 0.8 ± 2 

-4.3 ± 0.4 ± 1.1 

0.24 

-4.2 ± 4 ± 18 ± 7 

16.1 ± 2 ± 11 

11.0 ± 1 ± 7 

-7.3 ± 0.7 ± 4.5 

0.30 

-8.9 ± 4 ± 19 ± 4 

-26.6 ± 2 ± 14 

-22.8 ± 1.5 ± 11 

-18.8 ± 1.2 ± 9.3 

0.35 

9.6 ± 3 ± 12 ± 6 

-23.9 ± 2 ± 8 

-28.9 ± 2 ± 7.5 

-31.2 ± 1.5 ± 7.4 

0.42 

28.0 ± 5 ± 11 ± 7 

-13.9 ± 4 ± 9 

-26.6 ± 4 ± 10 

-37.9 ± 5 ± 12 

0.50 

36.3 ± 4 ± 17 ± 3 

-13.2 ± 4 ± 16 

-31.0 ± 5 ± 20 

-48.4 ± 6 ± 27 

0.60 

43.4 ± 3.5 ± 15 ± 4 

12.2 ± 3 ± 16 

-35.9 ± 4 ± 24 

-64.5 ± 7 ± 38 

0.70 

56.0 ± 3 ± 14 ± 6 

-0.1 ± 3 ± 18 

-28.4 ± 4 ± 30 

-71.7 ± 7 ± 53 

0.84 

69.0 ± 3 ± 13 ± 1.5 

15.3 ± 3 ± 19 

-8.7 ± 5 ± 36 

-48.4 ± 11 ± 74 

1.00 

85.3 ± 3 ± 11 ± 0.7 

25.7 ± 2.5 ± 17 

-7.0 ± 5 ± 37 

-81.1 ± 11 ± 84 

1.20 

94.2 ± 3.5 ± 10 ± 1 

53.7 ± 3 ± 17 

57 ± 7 ± 39 

62.5 ± 18 ± 101 

1.40 

102 ± 4 ± 11 ± 2 

68.6 ± 4 ± 20 

88 ± 7 ± 48 

123 ± 19 ± 133 

1.70 

114 ± 3 ± 16 ± 2 

92.9 ± 5 ± 20 

150 ± 11 ± 48 

295 ± 32 ± 142 

2.40 

120 ± 2.5 ± 9 ± 3 

108 ± 4 ± 16 

218 ± 14 ± 46 

572 ± 53 ± 152 

3.00 

124 ± 3 ± 8 ± 3 

107 ± 4 ± 10 



3.50 

113 ± 7 ± 18 ± 1 




4.20 

125 ± 4 ± 9 ± 3.5 

110 ± 4.5 ± 7 



5.00 

118 ± 5 ± 11 ± 4 

85.3 ± 7 ± 16 

153 ± 18 ± 59 

398 ± 61 ± 236 

6.00 

122 ± 5.5 ± 8 ± 2 

102 ± 6 ± 8 

219 ± 17 ± 18 

664 ± 84 ± 56 

8.40 


102 ± 4 ± 7 



10.00 

128 ± 11 ± 13 ± 4 



565 ± 85 ± 66 

15.50 

130 ± 3 ± 16 ± 4 

88.8 ± 3 ± 16 

187 ± 10 ± 30 

597 ± 51 ± 80 

30.00 

125 ± 4 ± 10 ± 2.5 

78.7 ± 5 ± 11 

158 ± 20 ± 23 



• finally, the deviations for each data set were 
summed in quadrature as independent values, 


W) = 


\ 


N s _ 

E( M n(i)(Q 2 )-K D (Q : 


( 26 ) 


where Ng is the number of available data sets. 
The resulting error is summed in quadrature with 
<5" orm (Q 2 ) to get the total systematic error on the 
n-th moment. 


The second type of error is related to the bin center¬ 
ing, interpolation and extrapolation. The bin centering 
systematic uncertainty was estimated as 

$n(Q 2 ) = y £ K n( x i,Q 2 )w i (Q 2 )5£ nt ( x i ,Q 2 ) , ( 27 ) 


where, according to the Nachtmann moment definition 
and the trapezoidal integration rule, one has 


K n (xi,Q 2 ) 

Wi(Q 2 ) 


: n+1 

^siObQ 2 ) 


Xj_ 

ii 


(*£2+1 1)/2 . 


n 2 M 2 x 2 & 

(n + 2) 2 Q 2 Xi 

( 28 ) 


The systematic error of the interpolation was esti¬ 
mated by considering the possible change of the fitting 


function slope in the interpolation interval, and was eval¬ 
uated as a difference in the normalization at different 
edges: 


<5s(Qo) 


\xi—x a \<A 

i 

|xj —Xfc|<A 

( 9i(xj,Q 2 0 ) 

3 3 


9i(xi,Ql)) 


9i ( x j 1 Q 0 )) 


( 29 ) 


where N and Nj are the number of points used to evalu¬ 
ate the sums. Since the structure function gi(x, Q 2 ) is a 
smooth function of x below resonances, on the limited re- 
interval (smaller than V) the linear approximation gives a 
good estimate. Thus, the error given in Eq. 11251) accounts 
for such a linear mismatch between the fitting function 
and the data on the interpolation interval. Meanwhile, 
the CLAS data cover all the resonance region and no in¬ 
terpolation was used there. The total systematic error 
introduced in the corresponding moment by the interpo¬ 
lation can therefore be estimated as 


/ cn+r 

dx —^-9f(x,Q 2 ) 

X a 

x n 2 M 2 x 2 £ 

£ (n + 2) 2 Q 2 x 


( 30 ) 















































IV. EXTRACTION OF LEADING AND 
HIGHER TWISTS 
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FIG. 6: Errors of the inelastic Nachtmann moment Mi: the 
empty circles represent statistical errors; the stars show the 
systematic error obtained in Eq. i HTTll : the low-x extrapolation 
error is indicated by filled squares. 


The systematic errors obtained by these procedures are 
then summed in quadrature to give, 


cw 2 ) = yjwmr + issmr + iw)) 2 • (si) 


In order to study the systematic error on the extrapo¬ 
lation at very low x we compared the moments extracted 
using different parameterizations of g\. We choose a 
Regge inspired form from Ref. Q and two QCD fits from 
Refs. Hal]- The difference was significant only for M \, 
for which the various errors are shown in Fig. Eland sep¬ 
arately given in Tabic ITTI 


According to Eq. m the contribution from the pro¬ 
ton elastic peak should be added to the inelastic mo¬ 
ments obtained above. The (^-dependence of the pro¬ 
ton elastic form factors is parameterized as in Ref. 26], 
modified accordingly to the recent data on Ge/Gm 37], 
as described in Ref. 38], The uncertainty on the form 
factors is taken to be equal to 3% according to the anal¬ 
ysis of Ref. 136}, and is added quadratically to both the 
statistic and the systematic errors. The elastic contribu¬ 
tion M®*(Q 2 ) turns out to be a quite small correction for 
Q 2 > n (GeV/c) 2 . Our final results for the total (inelas¬ 
tic + elastic) moments with n = 1,3,5 and 7 are shown 
in Fig. [3 Note also that the amount of the measured ex¬ 
perimental contribution to M n (Q 2 ) is at least 50%, and 
the systematic uncertainties increase significantly as Q 2 


In this section we present our analysis of the moments 
M n {Q 2 ) with n > 1. We extract both the leading and 
higher twist contributions to the moments, including a 
determination of the effective anomalous dimensions. 

Results for the first moment M\(Q 2 ) were presented in 
Ref. jnj. There the highest <5 2 -points [Q 2 > 5 (GeV/c) 2 ] 
were used to obtain the singlet axial charge, which 
for the renormalization group invariant definition in 
the MS scheme (which is adopted throughout this pa¬ 
per) gave: aQ nv = 0.145 ± 0.018(stat.) ± 0.103(sys.) ± 
0.041(low:r) ±o;dio (a s ), where the first and second errors 
are statistical and systematic, the third is from the x —> 0 
extrapolation, and the last is due to the uncertainty in 
a s . From the (^-dependence of the first moment the ma¬ 
trix elements of twist-4 operators were extracted, which 
allowed a precise determination of the color electric and 
magnetic polarizabilities of the proton (see Ref. 0 for 
details). 

As has been discussed in Refs. B1E10 , the ex¬ 
traction of higher twists at large x is sensitive to the 
effects of high-order pQCD corrections, for both the po¬ 
larized and unpolarized cases. In particular, the use of 
the next-to-leading order (NLO) approximation for the 
leading twist is known to lead to unreliable results for the 
determination of the higher twists in the proton F% at 
large x 19]. In this work we follow Refs, [j, 2,0], where 
the pQCD corrections beyond the NLO are estimated 
according to soft gluon resummation (SGR) techniques 
0 and a pure non-singlet (NS) evolution is assumed 
for n > 3 (H3- However, in contrast to Refs. [110, 
where SGR was considered for the quark coefficient func¬ 
tion only, we consistently add in this work the resumma¬ 
tion of large -71 logarithms appearing also in the one-loop 
and two-loop NS anomalous dimensions. This was previ¬ 
ously used in Ref. [ 20 ) to determine the strong coupling 
constant a s (M ^) from the experimental moments of the 
proton F -2 structure function determined in Ref. Q- 

Within the above framework, the Nachtmann moment 
of the leading twist part of the g 1 structure function, 
fign(Q 2 ), is (for n > 3) explicitly given by 


S Vn (Q 2 ) = 6A n [a s (Q 2 )f" S + e G ^ 

a s (Q 2 ) 


1 + 


47T 


( 


ocAnlo) . * ■ 

ZO DIS + ^7] 


l.NS)\ 
DIS ) 


(32) 


where the constant 8A n is defined to be the 71 -th moment 
of the leading twist at the renormalization scale / 1 2 , and 
7 ^ s is the one-loop NS anomalous dimension. In Eq. cm 
the quantity 5R^ S is given by 


SR 


NS 

n 


7/d(NLO) /-((NLO) ^(NLO) 
°DIS °n,LOG 

a (1,NS) _ a (1,NS) _ a (INS) 
^■In ^7dIS a 7„,LOG 


increases. 


(33) 
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FIG. 7: Total (inelastic + elastic) Nachtmann moments M n (Q 2 ) (filled circles) [see Eq. JT51 1 extracted from the proton world 
data in the range 0.17 < Q 2 < 30 (GeV/c) 2 for n = 1,3,5 and 7. Open squares and triangles correspond to the inelastic 
and elastic contributions, respectively. Statistical errors are reported for all the three terms; in case of the total moments the 
systematic errors are represented by the shaded bands. 


where 

A 7 f’ NS > E 7''’ NS > - (34) 

with 7 n* ;NS ’ ) being the two-loop NS anomalous dimension, 
/?o = 11 — 2Nf/3, (3i = 102 — 387V//3 and Nf the number 
of active quark flavors at the scale Q 2 . 

In Eq. 1331) 5C^ L0 ^ is the NLO part of the quark co¬ 
efficient function, which in the MS scheme is given by 


<5Cl NLO) = CpiS^n) 


S i(n) + | - -7-^—r 
2 n(n + 1 ) 


- &(n) + S + STT 


* - I 1 < 35 > 


where C F = ( N 2 - 1)/(2N C ) and S k (n) = Y!j=\ l/j k . 


For large n (corresponding to the large-a; region) the coef¬ 
ficient d NLO) is logarithmically divergent; indeed, since 
Si(n) = 7 E+log(n)+C , (l/n), where 7 e = 0.577216 is the 
Euler-Mascheroni constant, and S 2 (n) = 7r 2 /6 + 0(l/n), 
one gets 


= C^NLO) + c (NLO) + 0{l/n) ) (36) 


with 


^(NLO) 

°DIS 


Cf 


' 2 3 9 

If + 2 7e “ 2 



(37) 


and 


^(NLO) 

°n,LOG 


Cf ln(n) 


ln(n) + 27 E + - 


(38) 
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For the quantity Ayi 1 ’ 1 ^ in Eq. lETH) one obtains 

A 7 (1,NS) = A7 a,NS) + A7 (1,NS) + 0{l/n) 5 (39) 


where 


/y v ( 1 > NS ) 

^7dis 


and 


Cf 

Po 

C A 

Nf 

3^1 

Po 


A ^(!,NS) 

a 7n,LOG 


C F 

22 


27t 2 + 32S'(oo) — 45 i 3 (oo) — 


17 


-AV-msM- g 


47T 2 1 

”9” + 3 


+ 7 e 8 K - 4 


Cf 

Po 


SK - 4 W, 


ln(n) 


(40) 


(41) 


with Ca = N c , S( 00 ) = = 

-0.751286, 5 3 (oo) = 1.202057 and K = C A (67/18 - 
7t 2 /6) — 5./V//9. 

In Eq. (1321) the function G n (Q 2 ) is the key quantity 
of the soft gluon resummation. At next-to-leading log 
(NLL) accuracy one has 


Gn(Q 2 ) = In(n) Gr(A„) + G 2 (A„) + O 


ajln k ~\n) ,(42) 


where X n = Po a s (Q 2 ) ln(n)/4n and 

Gi(A) = Cf~f-t [A + (1 — A)ln(l — A)] , 


A.'y-fr + 3 8 j K 

G 2 (A) = -G F -^|±Aln(l-A)-G ° ln(l-A) 

PO Po 


C F ^ ln(l-A) 
Po 


1 + iln(l - A) 


(43) 


Note that the function G 2 (A) is divergent for A —» 1; 
this means that at large n ( i.e ., large x) SGR cannot be 
extended to arbitrarily low values of Q 2 . Therefore, to be 
sure that the SGR technique can be used reliably at NLL 
accuracy it is essential to check that A ra is small enough, 
which in our case means restricting the twist analysis to 
the Q 2 -range above 0.8 4- 1 (GeV/c) 2 . 

It is straightforward to see that in the limit A n <C 1 one 

has G n (Q 2 ) -» a s (Q 2 ) [2 G^^T At^Log]/ 4 ^ so that 
Eq. 13211 reduces to the well-known NLO approximation. 
This implies that adopting the usual two-loop approx¬ 
imation for the running coupling constant a s (Q 2 ), the 
twist -2 expression cut contains all the NLO effects and 
the resummation of all the large-n logarithms beyond the 
NLO. 

The different running of the leading twist induced by 
resummation effects beyond the NLO has been investi¬ 
gated in Ref. El for the unpolarized case, and in Ref. 0 


for the moments of the proton gi structure function. It 
was found that, with respect to the NLO approximation, 
SGR effects enhance significantly the Q 2 -evolution of the 
leading twist moments at Q 2 ~ few (GeV/c) 2 , and that 
such an enhancement increases as the order n of the mo¬ 
ment increases. 

As far as power corrections are concerned, several 
higher-twist operators exist and mix under the renor¬ 
malization group equations. Such mixings are rather 
involved and the number of mixing operators increases 
with the order n of the moment. A complete calcula¬ 
tion of the higher-twist anomalous dimensions is not yet 
available, and therefore one has to use specific models or 
some phenomenological ansatz. 

An interesting model for higher twists is the renor- 
malon model E, which can be used as a guide to esti¬ 
mate the a:-shape of the higher twists (or more precisely, 
of the twist-4 and twist -6 terms). The renormalon model 
contains only one free-parameter, which means that it 
predicts the dependence of the higher-twist contribution 
to the moments upon the order n up to an overall un¬ 
known constant. It is also characterized by the fact that 
the renormalon anomalous dimensions are the same as 
the leading twist ones. However, in Refs. Eiia h 
was already found that the renormalon model cannot ex¬ 
plain simultaneously the power corrections to the trans¬ 
verse and longitudinal channels. Moreover, several phe¬ 
nomenological extractions of hig her-twist anomalous di¬ 
mensions made in Refs. |7L IflL ll7L llfll l3^| suggest that the 
latter may differ significantly from the leading-twist ones. 
Therefore, in this work we use the same phenomenolog¬ 
ical ansatz as adopted in Refs. j7, 9) Il3, Djl E (and in 
Ref. E for the n = 1 moment), which does not exclude 
the renormalon picture, but is more general. 

To be specific, the Nachtmann moments are analyzed 
in terms of the following twist expansion: 

M^(Q 2 ) = 6 Vn (Q 2 ) + HT„(Q 2 ) , (44) 

where the higher-twist contribution HT n (Q 2 ) is com¬ 
prised of twist-4 and twist -6 terms of the form 


HT„(Q 2 ) 


+ 



■ g s (Q 2 ) 

_a s (g 2 ) \ Q 2 


~ a s (Q 2 ) 

_a s {g 2 ) 




(45) 


where the logarithmic pQCD evolution of the twist-K con¬ 
tribution is accounted for by the term [a s (Q 2 )] <57 "" > with 
an effective anomalous dimension 5"/^\ and the param- 
eter 6an represents the overall strength of the twist-K 
term at the renormalization scale /i 2 . 

In Eq. (1451) only twist-4 and twist -6 terms are included. 
In practice the number of higher-twist terms to be con¬ 
sidered is mainly governed by the Q-range of the anal¬ 
ysis. Indeed, as the latter is extended down to lower 
values of Q 2 , more higher-twist terms are expected to 
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contribute. Here we note that: i) the inclusion of twist- 
4 and twist-6 terms works well for Q 2 > 1 (GeV/c) 2 , 
as already found in the case of the unpolarized moments 
B Hi Il3 , and ii) our least-y 2 fitting procedure turns 
out to be sensitive to the presence of a twist-8 term only 
f° r Q 2 1$ 1 (GeV/c) 2 , where the resummation of high- 
order perturbative corrections may start to break down. 
Therefore, we limit ourselves to considering only twist-4 
and twist-6 terms in the analyses for Q 2 > 1 (GeV/c) 2 . 

All the unknown parameters, namely the twist-2 coef¬ 
ficient SA n , as well as the four higher-twist parameters 
8a n 4 "*, 8j ^, San' 1 and Sji 6 ^, are for each order n simul¬ 
taneously determined from a ^-minimization procedure 
in the Q 2 range between 1 and 30 (GeV/c) 2 . Changing 
the minimum Q 2 value down to 0.7 -f- 0.8 (GeV/c) 2 does 
not modify significantly the extracted values of the vari¬ 
ous twist parameters. On the other hand, increasing the 
minimum Q 2 up to 2 (GeV/c) 2 leads to quite large un¬ 
certainties in the values of the twist parameters, due to 
a large decrease in the number of data points. 

The strong coupling constant in this analysis has been 
chosen to be a s (M§) = 0.118, consistent with the twist 
analysis of the unpolarized moments made in Ref. §|. 
The (arbitrary) renormalization scale fi is set to g = 
1 GeV/c. We point out that the high-Q 2 subset of the 
unpolarized Nachtmann moments of Ref. [f!| were ana¬ 
lyzed in Ref. [§0j in order to extract the value of a s (M§), 
including SGR effects up to NLL accuracy. The value 
found, a s {M 2 ) = 0.1188 ± 0.0010(stat.) ± 0.0014(sys.) 
(or 0.1188±0.0017 adding the errors in quadrature), was 
in full agreement with the latest Particle Data Group 
world-average value a s (M§) = 0.1187 ± 0.0020 Eoj . 

The fitting procedure provides the best-fit values of 
the twist parameters together with their statistical uncer¬ 
tainties. The systematic uncertainties are, on the other 
hand, obtained by adding the systematic errors to the ex¬ 
perimental moments and repeating the twist extraction 
procedure. Our results, including the uncertainties for 
each twist term separately, are reported in Table ITTT1 and 
in Fig. 0 |H^. The ratio of the total higher-twist con¬ 
tribution, HT n (Q 2 ), to the leading twist term 8g n {Q 2 ), 
is shown in Fig. 0a). Note that since the leading twist 
component of the moments is directly extracted from the 
data, no specific functional shape for the leading twist 
parton distributions is assumed in our analysis. In the 
same way also our extracted higher twists do not rely 
upon any assumption about their a;-shape. 

Our main results for the higher twists in Figs. 00 can 
be summarized as follows: 

• The extracted twist-2 term yields an important 
contribution in the whole Q 2 -range of the present 
analysis; it is determined quite accurately with an 
uncertainty which does not exceed 15% (statistical) 
and 20% (systematic); 


pure 1/Q 2 dependence, or may even become nega¬ 
tive at the lowest values of Q 2 and large n. This 
requires in Eq. <1351) a twist-6 term with a sign op¬ 
posite to that of the twist-4. As already noted in 
Refs. 0E1II3, such opposite signs make the total 
higher-twist contribution smaller than its individ¬ 
ual terms (see dashed lines in Fig. 0 ; 

• The extracted values of the higher-twist anomalous 
dimensions appear to be significantly larger than 
the corresponding ones of the leading twist (viz., 
7 ^ s = 0.67,0.97,1.17 for n = 3,5,7, respectively, 
at Nf = 4); 

• The total higher-twist contribution is important 
for Q 2 « few (GeV/c) 2 , and is still non-negligible 
even at Q 2 ~ 10 (GeV/c) 2 for the higher moments. 
Comparison with the higher twists extracted from 
the moments of the unpolarized structure func¬ 
tion Q in Fig.0clearly shows that the total higher- 
twist contribution is significantly larger in the po¬ 
larized case, as already observed in Ref. 7) and also 
in agreement with the findings of Ref. [41j . 


The extracted twist-2 contribution is given in Table Hvl 
and in Fig. 03 where it is compared with several NLO 
parameterizations of s pin -dependent parton distribution 


For n = 1 the twist- 
agrees well at large Q 2 


functions (PDFs) 

2 moment obtained in Ref. 
with the results of Refs. |42|, [43[, whereas at lower Q 2 
our findings are below the predictions of all the four 
PDF sets. We should note, however, that in Ref. m 
a next-to-next-to-next-to-leading order (N 3 LO) approx¬ 
imation was adopted, since for the n = 1 moment the 
SGR effects are totally absent. This gives rise to a run¬ 
ning of the leading twist which is faster than that at 
NLO. As n increases, our extracted twist-2 runs faster 
around Q 2 ss few (GeV/c) 2 , in agreement with the 
findings of Refs, {ii, ;Tf|,i.e. the running is enhanced by 
SGR effects with respect to the NLO scheme adopted in 

Refs. HI EE El El. 


Note that at large Q 2 (> 10 (GeV/c) 2 ) the extracted 
twist-2 contributions for n > 1 in Fig.03is systematically 
below the parameterizations in Refs. [H, El El Ell j with 
the discrepancy increasing with the order n. This would 
imply PDFs lower than those of Refs. asms 
at large x. Such an effect may at least partially be due 
to the neglect, or a different treatment, of higher-twist 
effects in the analyses of Refs. (3 HE EE El . which 
were carried out in x-space (see e.g., Ref. EH)- To fully 
unravel the origin of the above differences is, however, 
beyond the aim of the present paper. 


V. CONCLUSIONS 


• The Q 2 -dependence of the data leaves room for a 
higher-twist contribution which runs slower than a 


We have presented a self-consistent analysis of world 
data on the proton g i structure function in the range 
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TABLE III: Leading twist Sr/n and higher-twist parameters, appearing in Eq. (Ii51l . extracted from the Nachtmann moments 
for n > 3 at the scale Q 2 = 1 (GeV/c) 2 . The first errors are statistical, while the upper and lower ones are systematic. 



m 3 

M s 

m 7 

Sg„ 

0.0147 ± 0.0005±o.oo23 

0.0057 ± 0.0008to.ooo7 

0.0038 ± 0.0005to.ooo2 

SaW 

0.020 ± O.OOlto.oo? 

0.0155 ± 0.0007to;ooo9 

0.0103 ± o.ooosjio'.Mie 


2.2 ± 0.3tg'j 

2.3±0.5l^ 

2.6 ± 0.4lg j 

] 

-0.012 ± 0.002t° °g® 

-0.0127 ± 0.0009±g;ggg| 

-0.0108 ± 0.0005±g;ggg| 


3.0 ± 0.6±?;| 

2.4±0.8j:^2 

2.9 ±0.51^2 




1 10 


Q 2 (GeV/c ) 2 




1 10 


Q 2 (GeV/c ) 2 


FIG. 8: Results of the twist analysis for n = 1 (adapted from Ref. (IlJ) and for n = 3, 5 and 7 obtained in this work. Open 
circles represent the Nachtmann moments, and the solid lines are fits to the moments using Eqs. mi . mi and mi with the 
parameters listed in Table ITUl The twist-2 (dotted), twist-4 (dot-dashed), twist-6 (triple-dot-dashed) and total higher twist 
(dashed) contributions are shown separately. The errors indicated are statistical. 


0.17 < Q 2 < 30 (GeV/c) 2 , including recent measure¬ 
ments performed with the CLAS detector at Jefferson 
Lab §|. This analysis has made it possible to accurately 
compute for the first time the low-order moments of g\ 
and study their evolution from small to large values of 


Q 2 . Our analysis includes the latest experimental results 
from Jefferson Lab for the ratio R = g^/gt and a new 
model for the transverse asymmetry Ai in the resonance 
production regions, as well as the unpolarized cross sec¬ 
tions measured recently in the resonance region at Jef- 
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TABLE IV: The extracted leading twist contribution g n (Q 2 ) [see Eq. 11531 1 . reported with statistical and systematic errors. 


Q 2 [(GeV/c) 2 ] 

6rn (QA 

<MQ") x io- a 

Sg 5 (Q“) x 10-" 

Sgr (Q") x 10““ 

1.00 

0.1127 ±0.0030 ±0.0109 

1.47 ±0.05 ±0.24 

0.57 ±0.08 ±0.08 

0.380 ± 0.052 ± 0.048 

1.20 

0.1148 ±0.0030 ±0.0109 

1.40 ±0.04 ±0.23 

0.47 ± 0.06 ± 0.06 

0.248 ±0.034 ±0.031 

1.40 

0.1162 ±0.0030 ±0.0108 

1.35 ±0.04 ±0.22 

0.41 ± 0.05 ± 0.05 

0.193 ±0.026 ±0.024 

1.70 

0.1176 ±0.0030 ±0.0108 

1.28 ±0.04 ±0.21 

0.36 ± 0.05 ± 0.05 

0.153 ±0.021 ±0.019 

2.40 

0.1195 ±0.0037 ±0.0135 

1.18 ±0.04 ±0.19 

0.29 ± 0.04 ± 0.04 

0.109 ±0.015 ±0.014 

3.00 

0.1203 ±0.0037 ±0.0134 

1.13 ±0.03 ±0.18 

0.27 ± 0.04 ± 0.04 

0.096 ±0.013 ±0.012 

3.50 

0.1208 ±0.0037 ±0.0134 

1.10 ±0.03 ±0.18 

0.25 ±0.03 ±0.03 

0.088 ±0.012 ±0.011 

4.20 

0.1213 ±0.0037 ±0.0133 

1.06 ±0.03 ±0.17 

0.24 ±0.03 ±0.03 

0.081 ±0.011 ±0.010 

5.00 

0.1217 ±0.0037 ±0.0133 

1.03 ±0.03 ±0.17 

0.23 ± 0.03 ± 0.03 

0.075 ±0.010 ±0.009 

6.00 

0.1222 ±0.0036 ±0.0133 

1.00 ±0.03 ±0.16 

0.21 ±0.03 ±0.03 

0.070 ±0.010 ±0.009 

8.40 

0.1229 ±0.0036 ±0.0132 

0.95 ±0.03 ±0.16 

0.20 ±0.03 ±0.03 

0.062 ±0.008 ±0.008 

10.00 

0.1232 ±0.0036 ±0.0132 

0.93 ±0.03 ±0.15 

0.19 ±0.02 ±0.02 

0.058 ±0.008 ±0.007 

15.50 

0.1239 ±0.0036 ±0.0132 

0.88 ±0.03 ±0.14 

0.17 ±0.02 ±0.02 

0.051 ±0.007 ±0.006 

30.00 

0.1247 ±0.0032 ±0.0115 

0.81 ±0.03 ±0.13 

0.15 ±0.02 ±0.02 

0.043 ± 0.006 ± 0.005 
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Within the framework of the operator product expan¬ 
sion, we have extracted from the experimental moments 
at Q 2 ^ 1 (GeV/c) 2 the contributions of both leading and 
higher twists. Effects from radiative corrections beyond 
the next-to-leading order have been taken into account 
by means of soft-gluon resummation techniques. 

The leading twist has been determined with good ac¬ 
curacy, allowing detailed comparisons to be made with 
various NLO polarized parton distribution functions ob¬ 
tained from global analyses in Bjorken-x space. A faster 
running in Q 2 is observed in our twist-2 moments due to 
the inclusion of resummation effects beyond NLO. The 
twist-2 moments are also found to lie slightly below those 
calculated from the standard polarized PDFs, suggesting 
that the latter overestimate the leading twist at large x. 
This may reflect the different treatment of higher-twist 
effects in our analysis compared with those in the global 
PDF fits. 

The contribution of higher twists to the polarized pro¬ 
ton structure function gi is found to be significantly 
larger than for the unpolarized proton structure function 
F 2 , although some cancellations between different twists 
occurs at low Q 2 . 

Improvements in the determination of both the lead¬ 
ing and higher twist terms are expected to come with 
the availability of new CL AS data taken at Jefferson Lab 
with the 6 GeV electron beam, which will provide an ex¬ 
tended kinematical coverage up to Q 2 « 5 (GeV/c) 2 . Be¬ 
yond this, we anticipate significant progress in the mea¬ 
surement of polarized structure functions at higher Q 2 
and over a larger range of x with the upgrade of the Jef¬ 
ferson Lab electron beam to 12 GeV. 
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APPENDIX A: FIT OF THE PROTON 
TRANSVERSE ASYMMETRY A 2 

The parameterization of A 2 is based on an estimate of 
the polarized transverse structure function gx by means 
of resonance-background separation, where the resonance 
part is taken from a constituent quark (CQ) model 53 . 
while the background is described according Wandzura- 
Wilczek (WW) prescription [if|. As normalization, we 
use the Burkhardt-Cottingham (BC) sum rule s 4§j, for 
each Q 2 value of the data. The BC sum rule implies that 

1 

dx g 2 (x, Q 2 ) = 0 (Al) 

for any Q 2 , where the integration includes also the elastic 
peak. 

In practice it is more convenient to work with the 
purely transverse structure function gx, which is defined 
as 


9 t(x,Q 2 ) = gi{x,Q 2 ) +g 2 {x,Q 2 ) (A2) 

Decomposing gx into leading twist, elastic and higher 
twist terms, we can write 

9t(x, Q 2 ) = g™ w (x,Q 2 ) + gx(Q 2 ) < 5(1 - x) 

+ 9t T ( x iQ 2 ) (^ 3 ) 


This work was supported by the Istituto Nazionale 
di Fisica Nucleare, the French Commissariat a l’Energie 
Atomique, French Centre National de la Recherche Sci- 


where the first term represents the (twist-2) WW relation 
(which is found to be a good approximation in DIS), the 
second term represents the elastic peak contribution, and 
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This ansatz is motivated partly by duality arguments 
as well as by recent findings in polarized structure func¬ 
tion studies, which suggest a picture in which the reso¬ 
nance peaks fluctuate around a smooth background ex¬ 
trapolated from the DIS regime. Clearly this model ne¬ 
glects the interference between resonances and the back¬ 
ground, which can play an important role in the total 
cross section. However, given the absence of experimental 
guidance (at least above the two-pion production thresh¬ 
old), this approach is the minimal one suitable for the 
present analysis. 

Using the WW relation 45], one can rewrite qr in 
Eq. m as 

9 t{x,Q 2 ) = [ —gi(y,Q 2 ) 

Jx y 

+ 9t{Q 2 ) <5(! - x) + g^ T (x,Q 2 ) .(A6) 

From the BC sum rule in Eq. m and the Fubini theo¬ 
rem S3 we then find: 

/ dx/g* T {x 1 Q 2 ) = g?{Q 2 )-g e l{Q 2 ) 

J X 

= 8M ?Uq> Gm ( Q2 Wm(Q 2 ) - G E (Q 2 )) (A7) 


where Ge{Q 2 ) and Gm(Q 2 ) are the Sachs proton electric 
and magnetic form factors. 

The WW term g^ w is calculated from the phenomeno¬ 
logical parameterization of gi given in Ref. [7|, which is 
known to work well also in the resonance region and at 
the photon point ( Q 2 = 0). Furthermore, target mass 
corrections are applied in order to remove the kinemati- 
cal effects of working at finite Q 2 , 



FIG. 9: (a) Ratio of the total higher-twist [see Eq. 1451 1 to 
the leading twist given in Eq. iEH. Dotted line - A'h (from 
Ref. [il l ); triple-dot-dashed line - M 3 ; dashed line - M 5 ; solid 
line - M 7 . (b) Ratio of the total higher twist to the leading 
twist obtained in the analysis of the unpolarized moments in 
Ref. 0. 


the third parameterizes the remaining (higher twist) part 
of <7 t- 

Next we make use of an ansatz which assumes that the 
first term in Eq. (IA3I) . gjf w (x,Q 2 ), is due to the back¬ 
ground contribution and the second term, g^p-(x,Q 2 ), 
contains only the resonance part of the total cross sec¬ 
tion, 


(A4) 

(A5) 


where r = \J 1 + AM 2 x 2 /Q 2 . The resonance part of gr 
is directly related to the longitudinal-transverse interfer¬ 
ence term of the resonance production cross section, 


g r T es (W,Q 2 ) = - j vM * a LT '{W,Q 2 ) (A9) 


47 ^ay^Q 2 


where 


LT‘ 


' (W, Q 2 ) = J2 


N* 


Wq* 


B(W) 


x S* 1/2 (Q 2 )A 1/2 (Q 2 ) . (A10) 


Here the sum runs over all nucleon excited states N*, 
B(W) is the unit-area resonance shape described in the 
relativistic Breit-Wigner approximation, 


B{W) 


WM res r res 


9t W (x,Q 2 ) = g^ 6 (x, Q 2 ) , 
9t T (x,Q 2 ) = gT S {x,Q 2 ) . 


7r 


(W2-M r 2 s )2 + M 2 s r2 


(All) 
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FIG. 10: The leading twist moments (open circles) extracted in the present analysis for n > 3 and in Ref._|lJ] for n = 1, 
compared with the corresponding moments of various parton distribution sets: dotted 0; triple-dot-dashed j35j: dashed {42|; 
solid 0. 


and q* is the 3-momentum transfer in the resonance rest 
frame, 


<? 


* 



W 2 — M 2 
4T 


Q 2 


1/2 


(A12) 


A 1 / 2 (Q 2 ) an< i £ 1 / 2 (Q 2 ) in CQ models depends strongly 
on the choice of the potential and other model parame¬ 
ters. In order to improve this description we apply the 
BC sum rule given in Eqs. m and (1X71) to the entire 
resonance part of g™ s - This amounts to modifying g™ s 
by multiplying it by a factor 


The helicity amplitude A 1 / 2 (Q 2 ) is relatively well 
known for the most prominent resonances, while the lon¬ 
gitudinal amplitude Si/ 2 (Q 2 ) is largely unexplored ex¬ 
perimentally, apart from the A(1232) resonance for which 
some data do exist. Theoretical predictions for these am¬ 
plitudes can be obtained from CQ models which success¬ 
fully describe resonance mass spectra and some trans¬ 
verse electromagnetic couplings. We use the CQ model 
from Ref. |3| for both the A 1 / 2 (Q 2 ) and Sx/^Q 2 ) am¬ 
plitudes in order to calculate g™ s in Eq. (El'. 

Unfortunately, the Q 2 -evolution of the couplings 


n( 0 2, = 9t\Q 2 )-9^Q 2 ) 

W J fg th dx g^ s (x, Q 2 ) ’ 


(A13) 


Therefore, at each given Q 2 the BC sum rule defines the 
total area of the resonance structure function g^ s . 

The asymmetry A 2 can then be directly related to gr 
according to 


A 2 (x,Q 2 ) 


y/Q 2 grjx^Q 2 ) 
v Fi(x, Q 2 ) 


(A14) 
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where Fi(x,Q 2 ) is the familiar unpolarized structure 
function. The final parameterization is shown in Fig. EH 
compared with calculations of the MAID model from 
Ref. . The MAID results represent a sum over a few 
exclusive channels which should be reliable when W is 
not very large. New experimental data on g 2 in the res¬ 
onance region at different Q 2 values are clearly needed. 



FIG. 11: Constituent quark model calculations of A 2 (W, Q 2 ) 
in comparison with the MAID model predictions |4fi| at Q 2 = 
1.3 (GeV/c) 2 : triangles show the calculations as described in 
the text; solid (n production), dashed (n and g) and dotted 
(n, r /, K A and A'E) lines represent MAID model calculations. 
The dot-dashed curve indicates the upper Softer limit on A 2 . 


( N ). Therefore, in the DIS kinematics, defined here as 
W > 2 GeV, the asymmetry A 2 can be estimated through 
the WW formula within the systematic uncertainty of 
S sys (A 2 ) = 4 x 10 -2 . However, taking into account target 
mass corrections, which affect the gr structure function 
also in the DIS region, one finally finds 5 sys (A 2 ) = 1.6 x 
10 -2 [see Fig. E]. 



FIG. 12: Weighted difference between the experimental A 2 
values and the WW prescription A^ 2 W including the target 
mass corrections. 


In the DIS region data from Refs. [22, '23,0, [sij sug¬ 
gest that A 2 is rather small, and can be described within 
the WW approach. In order to quantify the agreement 
and to estimate the systematic uncertainty, we plot in 
Fig.^Jthe weighted difference between the data and the 
WW prescription, 


APPENDIX B: KINEMATIC HIGHER TWISTS 

In order to estimate contribution of the kinematic 
twists appearing in the expansion of the CN moments, 
we extract from our data the inelastic part of the d 2 mo¬ 
ment, defined as 


AA 2 


*exp _ iffW 
£*2 A 2 . 


( A15 ) d 2 (Q 2 ) = J dx x 2 {3g T (x,Q 2 ) - gi(x,Q 2 )} , (Bl) 


where Sa 2 is the A^ xp statistical error. One sees that the 
mean value within errors is compatible with zero, and 
the error of 4 x 10~ 2 has been estimated according the 
formula 


S sys (A 2 ) = 


' N 


, $A 7 ( X iiQi) 


- 1/2 


cr AA 2 


(A16) 


where <ta a 2 is th e width of the A A 2 distribution and 
the sum runs over all available A 2 experimental points 


where the structure function gr{x,Q 2 ) is described in 
Appendix m The extracted values of d 2 (Q 2 ) are given 
in Table Ivl and shown in Fig. 1131 

The lowest twist component in d 2 is twist-3, although 
higher twists can also contribute to d 2 at low Q 2 . Note 
that only the inelastic part of d 2 is extracted; the elas¬ 
tic contribution has to be added separately for a twist 
analysis of d 2 . The results indicate that at high Q 2 the 
values of d 2 (Q 2 ) are consistent with a vanishing twist-3 
contribution. 




































TABLE V: The inelastic part of d 2 (Q 2 ) extracted from data 
(see text). The results are reported together with the statis¬ 
tical and systematic errors. 
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Q" [GeW] 

d 2 (Q^ ■ 10-° 

0.17 

3.7 ± 1.6 ± 2.1 

0.20 

3.9 ± 0.8 ± 2.5 

0.24 

4.9 ± 1 ± 3.6 

0.30 

8 ± 1 ± 4.8 

0.35 

9.3 ± 0.9 ± 5.4 

0.42 

10.2 ± 2 ± 6.3 

0.50 

12.3 ± 1.8 ± 8 

0.60 

14.4 ± 1.4 ± 9 

0.70 

14.6 ± 1.2 ± 9.6 

0.84 

14.4 ± 1.2 ± 10 

1.00 

14.4 ± 1 ± 11 

1.20 

11.6 ± 1.2 ± 11 

1.40 

10 ± 1.2 ± 11 

1.70 

6.8 ± 1.5 ± 11 

2.40 

3.7 ± 1.3 ± 12 

3.00 

2.9 ± 1 ± 12 

3.50 

3.9 ± 0.5 ± 17 

4.20 

1.4 ± 1.1 ± 13 

5.00 

3.5 ± 1.6 ± 15 

6.00 

1.3 ± 1.3 ± 15 

10.00 

1.7 ± 1.2 ± 18 

15.50 

0.6 ± 0.7 ± 22 

30.00 

0.3 ± 0.9 ± 30 
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